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We determine the conditions for the existence of non-transverse factorizing magnetic fields in 
general spin arrays with anisotropic XYZ couplings of arbitrary range. It is first shown that a 
uniform maximally aligned completely separable eigenstate can exist just for fields h s parallel to a 
principal plane and forming four straight lines in field space, with the alignment direction different 
from that of h s and determined by the anisotropy. Such state always becomes a non-degenerate 
ground state (GS) for sufficiently strong (yet finite) fields along these lines, in both ferromagnetic 
(FM) and antiferromagnetic (AFM) type systems. In AFM chains, this field coexists with the non- 
transverse factorizing field h' s associated with a degenerate Neel-type separable GS, which is shown 
to arise at a level crossing in a finite chain. It is also demonstrated for arbitrary spin that pairwise 
entanglement reaches full range in the vicinity of both h s and h' s , vanishing at h s but approaching 
small yet finite side-limits at h' s , which are analytically determined. The behavior of the block 
entropy and entanglement spectrum in their vicinity is also analyzed. 
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I. INTRODUCTION 


The ground state (GS) of strongly interacting spin sys¬ 
tems immersed in a magnetic field h can exhibit, under 
certain conditions, the remarkable phenomenon of fac¬ 
torization ,H, i.e., of becoming a product of single spin 
states. Such exact factorization can occur at finite fields 
despite the strong couplings existing between the spins, 
albeit at very specific values (and orientation) of the field. 
In the seminal work of ref. Tj, it was shown that anti¬ 
ferromagnetic (AFM) chains with first neighbor XYZ 
couplings possess a separable Neel-type GS (NGS) if the 
field vector lies on the surface of an ellipsoid determined 
by the couplings. Factorization was then investigated in 
other models with transverse fields j2l-l20t , with a general 
formalism for describing factorization introduced and dis¬ 
cussed in 



we have shown that in finite XYZ 
chains, the transverse factorizing field (TFF) h zs point¬ 
ing along a principal axis ( z ) corresponds actually to 
the last GS S z parity transition (level crossing). The 
ensuing separable GS is two-folcl degenerate, breaking 
a basic symmetry of the Hamiltonian ( S z parity). The 
non-transverse factorizing fields (NTFF) h' s of [l| will be 
shown to also belong to this class in finite cyclic chains, 
i.e., they arise at the last GS level crossing and determine 
a degenerate separable GS breaking translational invari¬ 
ance (TI). In finite systems the underlying mechanism of 
factorization in these cases is the existence of separable 
linear combinations of the symmetry preserving entan¬ 
gled crossing states. 

In this work we first determine the general conditions 
for exact factorization under non-transverse fields. It is 
then shown that a uniform non-degenerate separable GS 
(UGS) does exist at a field h s which does not belong 
in general to the ellipsoid of ref. [l|, and does not cor¬ 
respond to a level crossing. This GS actually arises in 
both AFM and FM-type systems, even for couplings of 
arbitrary range provided there is a fixed anisotropy ratio, 


but only for fields parallel to a principal plane, with the 
set of fields h s forming four straight lines. Factorization 
emerges here from the splitting of the degenerate sepa¬ 
rable eigenstates existing at the TFF h zs . Unlike h ' s , h s 
can be arbitrarily strong, allowing the separation of the 
UGS from the remaining spectrum. This enables an easy 
preparation of an exactly separable state, which can be 
useful for quantum information applications (a product 
initial state is assumed in the standard model of quantum 
computation [2lj). 

A second but not less important aspect of factoriza¬ 
tion is that it corresponds to an entanglement transi¬ 
tion: In the transverse case, the factorizing field is, re¬ 
markably, the point where pairwise entanglement reaches 
full range in its immediate vicinity, and changes its type 
1 1 Hi ED, 0- We had previously shown that the 
entanglement between any two spins reaches in a finite 
chain weak yet finite common side-limits at the trans¬ 
verse field h zs , irrespective of separation or coupling 
range 0U3, arising from the entangled crossing states. 
This type of limit also occurs at the NTFF h' s of ref. [l[ , 
as will be shown. But in addition, we will prove that pair¬ 
wise entanglement also reaches full range at the vicinity 
of the NTFF h s leading to a non-degenerate UGS. Here 
the entanglement between any two spins, though 0 at h s , 
is turned on as h s is approached from either side, with 
the concurrence vanishing then linearly with | h — h s |. 
The underlying reason is essentially the monogamy of 
entanglement [22l . [23j | , which prevents distant pairs from 
becoming entangled if first or close neighbors are strongly 
entangled. In the vicinity of h s , close neighbor entangle¬ 
ment decreases strongly, allowing the emergence of weak 
yet non-zero entanglement between distant pairs. The 
behavior of the block entanglement entropy in the vicin¬ 
ity of the NTFF will be also analyzed. It will be shown to 
vanish essentially quadratically at h s , while at h' s it will 
approach finite side-limits in a finite chain, which will 
be analytically determined. The entanglement spectrum 
will indicate, as expected, just one nonzero eigenvalue at 
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h Sl although at the side-limits of h' s two nonzero eigen¬ 
values will remain. 

The general equations for non-transverse factorizing 
fields and its uniform and Neel-type solutions are de¬ 
rived and discussed in section El whereas entanglement 
together with illustrative results for the pairwise concur¬ 
rence, block entropy, entanglement spectrum and mag¬ 
netization in FM and AFM chains with XY and XYZ 
couplings under non-transverse fields are discussed in sec¬ 
tion III. Conclusions are derived in IV. 


II. FACTORIZATION IN NON-TRANSVERSE 
FIELDS 

A. General Equations 

We consider an array of n spins Si not necessarily 
equal, interacting through XYZ Heisenberg couplings of 
arbitrary range in the presence of a general magnetic field 
h l = (h x: hjhl), not necessarily uniform. The Hamilto¬ 
nian reads 

H = ~Y,K s i~\ E ■Wsj*, (!) 

where i,j label the sites in the array, S9, p, = x, y, z , the 
spin components at site i and Jj the coupling strengths 
between spins i and j (Jj > 0 corresponds to the FM 
case whilst J'j < 0 to the AFM case). In the transverse 
case h l x = h l y = 0 V i, H conserves the Sb-parity P z = 
exp[*7r ^2i(Si + &)] ([H,P Z ] = 0). This symmetry no 
longer holds for non-transverse fields. 

We now determine the general conditions for which H 
possesses a completely separable eigenstate of the form 

|0) = ®" = i-Ri|0i), Ri = exp [-ifaSf] exp [~i0iS ?], (2) 

where |0j) denotes the local state with maximum spin 
along z (S?|0i) = Sj|0i)) and Ri rotates this state to di¬ 
rection rii = (sindi cos fa, sin 9i sin</>j, cos 0i). The equa¬ 
tion H |0) = Eel©) leads, after writing H in terms of 
the rotated spins S9 = R^S'- Ft]. to: 

I) The field independent equations 

jj (cos fa cos fa — cos 0i sin fa cos 9j sin fa) = 

J l j (cos 6i cos (j>i cos 9j cos (fj — sin fa sin fa) 

+jj sin 9i sin 9j , (3) 

Jy (cos 9i sin fa cos fa + cos fa cos 9j sin fa ) = 

J l j (cos 9i cos fa sin fa + sin fa cos 9j cos fa ), (4) 

which are also independent of spin and are responsible 
for cancelling all elements connecting |0) with two-spin 
excitations, and 


II) The field dependent equations 

hi sin 9i — cos 9i(h l x cos fa + h l y sin fa) = 

Y2 Sj [ cos Si sin 9j (J l f cos fa cos fa + Jlj sin fa sin fa ) 

— J l j sin 9i cos 9j ], (5) 

h l x sin fa — hy cos fa = 

Y2 Sj sin 9j \—J x sin fa cos fa + Jj cos fa sin fa ], (6) 

which cancel all elements connecting |0) with single spin 
excitations and are just the mean field stationary equa¬ 
tions dgfaH) = 0, d ( p i (H) = 0, where 

(H) = (&\H\e) = - Y,(Si) ■ + 5 E Jij < S i ))’ W 

* j 

with (St) = Sjrij and J 13 a diagonal matrix of elements 
Jj. If Eqs. © © are satisfied V i.j. Eqs. © © de¬ 
termine the set of factorizing fields. 

In terms of the alignment directions n, = nf and the 

orthogonal unit vectors nf = (— sin fa, cos fa, 0), nf = 
(cos 9i cos (fi , cos 9i sin fa, — sin 9j ), we may also express 
Eqs. © © as 

nf ■ J'-'nj = nf ■ J"n'f , (8) 

nf ■ J' J nf = nf ■ J'-'nf , (9) 

which imply = jj, y , and J l j y , = -Jj, x , when writ¬ 
ing the coupling in © in terms of the rotated spins Sf , 
i e„ And Eqs. © © 

become 

nf ■(h i +^2J ij (S j )) = 0, p'=x',y', (10) 

3 

implying that h 1 should cancel the components of 
Y2j J 1 '' (Sj) orthogonal to the alignment direction, such 
that 

h 1 +Jfj ij (Sj) ocni. 

3 

The general solution for the NTFF at site i is then 

K = h \| + hf , ( 11 ) 

where h]^ = h ]| n^ is an arbitrary field parallel to the local 
alignment direction, which just shifts the energy 0, and 

hi = - ^2[J ij (Sj) - mim ■ J ij (Sj ))], ( 12 ) 

3 

is a field orthogonal to the alignment direction (n; • 
h‘ ± = 0), representing the NTFF of lowest magnitude. 
Nonetheless, a finite hi will be normally required in or¬ 
der that |0) be a GS (see sec. II.C). Let us remark, fi¬ 
nally, that Eqs. © m remain valid for general cou¬ 
plings E^W/inQ. 
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FIG. 1. (Color online) Schematic plot of the uniform solution. 
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B. Uniform solution 


Eqs. (0 | 6 ]) (or © m) are quite general and de¬ 
scribe a wide range of interesting scenarios. We examine 
first the possibility of a uniform solution with Oi = 0, 
(f>i = 0 V * (Fig. ©, such that |0) is a maximum spin 
state: |(]TT S))| = JA Si- Such solution preserves TI and 
then has the possibility to be a non-degenerate GS in 
systems with this invariance under a uniform field. 

Eq. © becomes (Jd — Jd ) cos 0 sin 20 = 0, implying, 
if Jd — Jd / 0 for at least one pair, that the spin vector 
(Si) should be parallel to a principal plane ( xz if 0 = 0 , 
yz if <f = 7 r /2 and xy if 0 = 7 r/ 2 ). Without loss of 
generality, we can assume 0 = 0 (the other choices are 
rotations of this case). Eq. © then leads to 


cos 2 0 


jij _ jij 

J y ° z 
jij jij 
z 


= x, 


(13) 


FIG. 2. (Color online) Factorizing fields for ferromagnetic 
(left) and antiferromagnetic (right) XYZ chains in the xz 
principal plane of the field space. Solid straight lines depict 
the fields determining a uniform ground state (UGS), whereas 
the dashed straight lines those determining uniform excited 
eigenstates (UES). The ellipse depicts the fields correspond¬ 
ing to the Neel-type ground state (NGS, solid lines) or ex¬ 
cited eigenstate (NES, dashed line). The plot corresponds 
to J 2 = 0 and J x > 0 (< 0) in the FM (AFM) case, with 
0 < J y / Jx < 1. The arrow indicates a direction of the exter¬ 
nal field along which one (FM) or two (AFM) GS factorizing 
fields are encountered as its magnitude increases. The field 
direction n 1 differs from the spin alignment direction ng . The 
insets depict the decomposition m of the non-transverse fac¬ 
torizing field for the UGS in both diagrams, with the dashed 
arrow indicating the transverse factorizing field h zs . 


For 0 = 0, Eqs. © (or (fTUl) ) imply h l y = 0, i.e., the 
field at each site should be parallel to the corresponding 
principal plane (xz). Eq. © then becomes 


if Jd zfz JJ?, implying a constant anisotropy ratio x for 
these pairs, and an isotropic coupling Jd = J n V /i if 
,/* J = ,/d. The condition 0 < y < 1 imposes the restric¬ 
tion 


rj > J l y j > d'j or Jd < Jd < Jd . ( 14 ) 

Eqs. c© © entail that the Jd should be of the form 
jd = jd +r dj /i) (15 ) 

with the s satisfying d. The state |@) will then 
depend just on y = J , v _ , being independent of the cou- 

Jx J z 

pling range determined by J lJ and r lJ . Notice that Eq. 
m leads to four possible alignment directions in the xz 
plane, corresponding to the solutions ±9 and ±( 7 r — 9 ), 
with 9 e (0, 7 t/ 2 ). 

We remark that in the fully isotropic case nj = 0 V j 
in (1151) (rotationally invariant coupling), 9 and 0 remain 
obviously arbitrary under Eqs. © ©, whereas in the 
XX case Jd = Jd V ij (coupling invariant under any 
rotation around the z axis), Eq. a is trivially satisfied 
while © leads to sin0 = 0 if Jd zfz J l _j for at least one 
pair, in agreement with m, implying alignment just 
in the z direction. We will focus in what follows on the 
anisotropic case 0 < y < 1 , where the alignment direction 
is non-trivial (9 £ ( 0 , 7 r/ 2 )). 


h l z sin 9 — h l x cos 9 = h\ , (16) 

where 

hf = sin 9 cos 9 $j (4’ - 40 • (3.7) 

j^i 

Setting ng = rii = (sin0,0,cos0) and n0 = 
(—cos 0,0, sin 0), Eqs. (flGl) (fT71) imply that the NTFF 
is given by 

h l s = h^ng + h\rig , (18) 

in agreement with m , with h|| arbitrary and h\njf or¬ 
thogonal to the alignment direction. Eqs. d © give 
rise to a family of NTFF lying along four straight lines 
(Fig.©, one for each alignment direction. 

Note that field and spin directions cannot be parallel 
if h l ± ^ 0: At fixed field direction n 7 = (sin 7 , 0, cosy), 
i.e., h\ = h l s (y)n 7 , Eqs. (fT51) leads to 


h i 

Kii) = 


sin (0 — 7 ) 


(19) 


which diverges for 7 -j 9. When the four values of 9 are 
considered, Eq. (fTHl) leads to two distinct values of \h l s \ 
at fixed 7 7 ^ ±0, which merge at the principal axes (Fig. 

©. 
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For 7 = 0 , we recover from (HU) the TFF 0, 0 


Ks = K(o) = 


sin 6 


( 20 ) 


which is the solution of (fTBl) for h x = 0. We can then 
also express Eq. (fl8l) as (n z = (0,0,1)) 

K = h l n g + h\ s n z , (21) 


where h 1 = hi — h\ / tan 9. Hence, we can also consider h\ 
as the sum of the TFF h\ s = h l zs n z plus a non-transverse 
field of arbitrary magnitude hi along the spin alignment 
direction n g , which just shifts the energy Eq. 

In systems with TI (i.e., infinite or cyclic), Si = S and 
h l _|_ = h_ l V i, implying a uniform factorizing field h s ( 7 ) 
at fixed orientation 7. Nonetheless, Eqs. mi) cu) show 
that the uniform solution remains feasible even in the 
absence of TI, provided the h 1 at each site can be con¬ 
trolled independently. In particular, in open finite uni¬ 
form chains or lattices with short range couplings, the 
uniform separable solution requires just border correc¬ 
tions to the otherwise uniform bulk factorizing field. 


H if hy = 0 and h l x > 0 V i. A GS with expansion coef¬ 
ficients real and of the same sign in this basis will then 
exist, as different signs will not decrease ( H). But such 
GS cannot be orthogonal to |0) if 9 £ ( 0 , 7 r) (implying 
h l > 0 in (l 2 lT) if h x > 0 ), so that it must coincide with 
|0) when |0) is an exact eigenstate. The case h l x < 0 can 
be reduced to the previous one by a rotation of angle n 
around the z axis, which leaves the rest of H unchanged. 

Besides, in the transverse case hi = 0 Vi, the states 
| 0 ) and | — 0 ) = -fk| 0 ), obtained for 6 = ±| 0 |, become 
degenerate (Eq. E51) i. The TFF hl s determines then a 
pair of degenerate UGS | ± 0) when EH) holds 0 , and 
the addition of a field parallel to ng (n_g) removes this 
degeneracy, leaving just |0) (| —0)) as GS. The transition 
to the UGS takes then place at h l zs . □ 

The gap to the first excited state can then be made 
arbitrarily large by increasing the fields hi (Eq. E3l) h 
Note that the similar case > | Iff \ V i,j can be reduced 
to the previous one after a 7r/2 rotation around the y 
axis. Hence, in this case the transition takes place at the 
transverse field along x, h xs = h l ± / cos 9 = h l zs tan 0 . 


D. Neel-type solutions. 


C. Uniform ground state 

For the uniform solution, the energy 0 becomes 

Eo=- W 4 7 ' ~ J y j + J l j ) - E S ^\\ ( 22 ) 

hj i 

= -!E + J y J i j ) E s ' h ' ■ (23) 

i,j i 

It is then apparent that |0) will be GS if the fields h'^ng 
(or equivalently h l ng) along the spin alignment direc¬ 
tion are sufficiently strong , since no other state has an 
energy which decreases more rapidly with the applied 
field. Therefore, a transition to this uniform separable 
GS (UGS) will always arise as hi increases, in both FM 
or AFM-type systems, as can be appreciated in Fig. [ 2 ] 
(transition from dashed to solid along the straight lines). 
Before this transition, |0) is an excited eigenstate (no 
other state can increase its energy more rapidly for de¬ 
creasing h l, s). 

We now show that if Eq. m is satisfied and, V i,j, 

4 j > \4 j \ > ( 24 ) 

such transition occurs at the TFF h l zs , i.e., |0) will be 
GS V hi > 0 in ET1) (left panel in Fig. [ 2 ]). 

Proof: We first note that if <f>i = 0 and 9i = 6 V i, Eq. 
© leads to 

|0) = ® i (E^o( 2 f) 1/ “ c °s 2S *- fc fs in El^))> where 

Sf\ki) = (Si — k)\ki). Eq. Ell) implies that the interac¬ 
tion in FI will contain just negative or zero off-diagonal 
elements in the standard basis (®i|fci)}, as seen by writ¬ 
ing JTJ in terms of Sf = Sf ± iSf. The same holds for 


In addition to the uniform solution, other solutions 
of Eqs. (0 0 can exist, which break TI. This is the 
case of the Neel-type separable eigenstates determined 
in PJ for the AFM chain with first neighbor couplings in 
a uniform field (Jy = 0 , = S iy j±i in (fTH) . with J^ < 0 

for /i = x,y,z ), where 0i,(f>i have alternating values. In 
a finite cyclic chain (with an even number n of spins), 
such solution must then be two-fold degenerate, arising 
at the crossing of two non-separable TI eigenstates. The 
mechanism is then similar to that of the TFF for the 
uniform solution 0. The associated NTFF h' s points 
to the surface of an ellipsoid jl|, given for S = 1/2 by 

Kl < 

(Jr + Jy) (Jr + Jz) (Jy + Jz)(Jy + Jr) 

h' 2 

__ = 1 Cot;'! 

(Jz + Jx)(Jz + Jy) 1 

Within the xz plane, h' s = /i(.( 7 )n 7 describes an ellipse 
(Fig. 0, satisfying 

r — ■ \ 12 (Jr + Jz)(Jx + Jy)(Jz + Jy) / 0R \ 

(Jr + Jy) COS 1 7 + (Jz + Jy) Sin 7 

While in a FM-type chain such solution also exists but 
corresponds to an excited eigenstate (left panel in Fig. 
0 . in the AFM case it is a GS which coexists with the 
previous UGS in the xz field plane (right panel). For 
instance, they can arise for the same field orientation 
at different field magnitudes. This possibility is related 
with the existence of different solutions for the local uni¬ 
tary operations which can leave an eigenstate invariant 
in the treatment of M- Moreover, the point where the 
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FIG. 3. (Color online) Factorization diagram (top panels) 
and factorizing fields (central and bottom panels) for XYZ 
couplings satisfying Eq. m, with J x > 0 (< 0) in the left 
(right) panels and J y of the same sign as J x . The different 
combinations of couplings are indicated (see also text). For 
Jx < J y < 0, the Neel separable eigenstate ceases to be GS 
when J z > —J y , although the uniform separable eigenstate 
remains GS for appropriate fields, as seen in the top and bot¬ 
tom right panels. Eq. (12511 may determine an hyperboloid 
when the couplings have different signs, as seen in the top 
and bottom panels. 


straight line of the uniform solution crosses the ellipsoid 
( h s = h' s ) is precisely that beyond which the uniform 
solution becomes GS (at this point the Neel-type solu¬ 
tion becomes uniform, coinciding with the UGS). Hence, 
within the first quadrant, the UGS arises for field angles 
0 < 7 < 9 in the FM case but 9 < 7 < 7 m in the AFM 
case, with 


tanq m 


j T Jy 
Jy T Jz 


tan 9. 


(27) 


of couplings for the case m are considered. When \J Z \ 
increases from 0, the diagrams of Fig.[2]remain essentially 
unchanged if \J Z \ < \J y \ (central panels in Fig. [3]), both 
in the proper FM and AFM cases (all couplings of the 
same sign) as well as in those where J z has the opposite 
sign of J x . However, when \J Z \ > \J y \ (with (ITU) still 
holding, e.g. J x < J v < 0 and J z > —Jy), the ellipsoid 
turns into an hyperboloid and the Neel-type state ceases 
to be GS in the originally AFM case ( J x and J y negative), 
as indicated in the top and bottom right panels. Yet the 
uniform separable eigenstate remains GS for both J x > 0 
and J x < 0 (blue lines in the bottom panels). This is still 
the case when | J z | increases beyond | J x |, as indicated in 
the top panels. We just mention that the cases J x > 
0 > J y > J z and J x < 0 < J y < J z are equivalent, 
respectively, to J z > 0 > J v > J x and J x > J y > 0 > 
J z , after rotation around the y axis. Furthermore, cases 
where Eq. ED does not hold can be transformed to the 
present situation by a suitable rotation. 


III. ENTANGLEMENT IN THE VICINITY OF 
FACTORIZATION 


A. Entanglement in the vicinity of the UGS 


Let us now discuss entanglement in the vicinity of the 
NTFF h s leading to the uniform GS |0). For simplic¬ 
ity we consider here a uniform field h in a spin S sys¬ 
tem with TI, where the reduced two-spin density matrix 
Pij = Trjy|GS')(GS'| (Trjy denotes the trace over the 
complementary subsystem) depends just on the separa¬ 
tion between the two spins. This reduced state will be in 
general a mixed state when | GS) is entangled. And such 
mixed state is said to be entangled if it cannot be written 
as a convex mixture of product states pi ® pj [IH, i- e -, if 
it cannot be generated by local operations and classical 
communication [ 2 l| . 

We first show that pairwise entanglement reaches full 
range in the vicinity of the factorizing field h s . 

Proof. For h close to h s , the GS can be obtained by 
considering first order perturbative corrections to |@): 


\gs) « |9) + •(E < g i )|e) | t/) 

„ E v — Eq 


— | 0 ) + (aX)5, Sj + .. .)|0)(28) 

i i,j 


Within this window, the GS of the AFM chain exhibits 
then two distinct factorizing fields as h increases at fixed 
n 7 (right panel in Fig. [2D, a result which has not yet been 
reported. 

Eq. (125[) may also determine an hyperboloid when J x , 
J y and J z do not have all the same signs, as shown in 
Fig. [3] where all (non-equivalent) possible combinations 


where \v) are the exact excited eigenstates at h s (H\v) = 

_ / 

E v \v), (i/|0) = 0), normally entangled, and S i = 
RiS~R\ the rotated lowering operators, with a, faj 
and all remaining terms of order Sh± if h — h s = 
5h±njp + 6h\\ng. In the rotated standard basis 
(S l z ,\k'i) = (Si — k)\k'i)) and considering first S = 1/2, 
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Eq. (l28l) leads to 


Pij 


f 1 

a 

a 

Pi. 

a 

0 

0 

0 

a 

0 

0 

0 

KPij 

0 

0 

0 


+ 0(6h 2 ± ). 


just one nonzero eigenvalue p\ = 1 , whereas in its vicinity 
the remaining eigenvalues will be small, of order 0(<5/i^). 


(29) 


B. Entanglement in the vicinity of the NGS 


According to the positive partial transposition criterion 
[25l [ 2 (| , this state will be entangled if its partial trans¬ 
pose pF is non-positive, i.e., if it has at least one neg¬ 
ative eigenvalue. But the the partial transpose of m 
has eigenvalues 1,0 and ±/3y up to 0(Sh±), so that py 
will be entangled if by 7 ^ 0. And the exact coefficients 
{3ij obtained from (12811 are not strictly zero for any pair 
i. j linked by successive applications of the couplings in 
H , due to the two spin excitations present in the exact 
eigenstates \v). 

For higher spins S, pij will be more complex (of 
( 2 S + l ) 2 x ( 2 S + l) 2 ) but will still contain a first sub¬ 
matrix of the form H29\) . Hence, it will also be entangled 

if pij 7 ^ 0 , since the partial transpose of this block is 

T 

the first block of the full partial transpose pF, and is 
non-positive at 0(6h±). This prevents the full pF from 
being positive semidefinite (in which case all principal 
submatrices should also be so). □ 

For S = 1/2, the entanglement between spins i and 
j can be measured through the concurrence [27] Cy = 
2A max — TV My, where A max is the largest eigenvalue 

of the matrix My = [p]- 2 pijp]- 2 ] 1 / 2 , with ,5y = cr y <g) 
ijyP*jCry (g) Uy in the standard basis. Up to 0(6hj_) , Eq. 
m then leads to 


Cij k, 21 Pij | oc \Sh.±\ 


(30) 


Note that at this order, a in o has no effect on the 
eigenvalues of pj- nor on Cij. Eq. (l3()l) implies that Cy, 
while acquiring finite positive values in the neighborhood 
of h s , will vanish linearly (as \8h±\) as ft —> h s , i.e., as 
it crosses the straight line of factorizing fields at a fixed 
direction n 7 . The corresponding entanglement of forma- 

tion [27(|, Eij = - lo §2 Pu, withp ± = — v 2 " ij , 

is just a convex increasing function of Cy, which vanishes 
as — A Cfj log 2 (C 2 /4e) for Cy —> 0. Hence, for ft —» h s it 
will vanish essentially as —8h\ log 2 \8h j_|. 

It is also seen from (EH1) that the eigenvalues of py will 
be either 1 (with negative 0(8hf\_) corrections) or small 
( 0 ( 6 / 1 ^)). Hence, the entropy S'(py) = —Trpy log 2 Py , 
which measures the entanglement between the pair and 
the rest of the system [28], will also vanish essentially 
as — 5h\ log 2 |<5ftj_| for ft —> h s . The same behavior at 
h s will be exhibited by the single spin entropy S[p(l)], 
where p(l) = pi = Trjpy denotes the single spin reduced 
state, and also by the block entropy [29| S[p(m)] of m 
contiguous spins, where p(m) denotes their reduced state. 
Factorization can in fact be directly seen through the 
entanglement spectrum [ijl,i3C|, i.e. the set of eigenvalues 
of the reduced states p(m). At h = h s . p{m) will have 


We first recall that in the transverse case h = hn Zl 
the behavior of Cy close to the TFF h zs in the GS of 
a finite FM-type chain mm is different from that de¬ 
scribed above. Since in the transverse case the S z par¬ 
ity P z is conserved, the exact GS of a finite spin chain 
has a definite parity, exhibiting parity transitions (the 
last one at the TFF h zs ) as the transverse field is in¬ 
creased [lot . This implies that for h —> hf s , it actually 
approaches the entangled definite parity degenerate side- 

limits |© =l= ) = ,l e > =1 =l~ e > with I — 0) = P,|0). These 

1 1 - v /2(i=t<—e|e))’ 1 1 z ' ' 

states lead to common finite side-limits C± of the con¬ 
currence Cy for any pair i 7 ^ j, given for S = 1/2 by 

m 


± sin 2 9 cos” 2 9 
1 ± cos ”9 


(31) 


where n is the number of spins and cos 9 is determined 
by UM, with (—0|0) = cos” 9. For finite n, a small 
but finite discontinuity in Cy is then encountered as the 
transverse field h crosses h zs , reflecting the parity change 
of the GS at h zs . Of course, exactly at h = h zs , the 
GS is two-fold degenerate and entanglement depends on 
the choice of GS, as in general degenerate systems f3l| . 
Factorization implies that the minimum entanglement at 
this point is zero (obtained when choosing | ± 0) as GS), 
even though the side-limits are finite. 

Remarkably, in the AFM chain, Eq. m remains for¬ 
mally valid for the side-limits of Cij at the Neel NTFF 
h' s , i.e., as ft. at a fixed orientation n 7 crosses the ellipsoid 
of factorizing fields h' s . The reason is that the exact GS 
of a finite cyclic AFM chain in a uniform field preserves 
TI away from crossing points and hence, it approaches 
for ft — > h'f the entangled TI side-limits 


i0±\ _ |0jy) ± | - Qjy) 

~ N \/2(l ± (—0 jv|0jv)) 


(32) 


where \<d N ) = | 0 i</>i, 9 2 p2, ■ ■ •), | - ©at) = 

| 02 </ 2 ,...) = T\Q n ) denote the degenerate 

Neel-type separable GS’s at h' s (T denotes the one-site 
translation operator, with T|0^) = ±|0^)). And 
these states | 0 ^j) lead to similar side-limits for the 
concurrence Cy between any two spins i j (see (|35l) 
below), i.e., 


q ± _ I s i n2 O' cosn 2 O' 
1 ± cos" 9' 


(33) 


where 9' is half the difference between the alternating 
angles of the Neel solution. This angle is determined by 
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a 

2 f\l ( Jz~\~Jy)(Jx~\~Jy ) ( Jx~\~Jy ) COS 2 4~ Jy) sin 2 7 

<<)h Jx+Jz {Jx + Jy) 2 COS 2 J+iJz+Jy) 2 sin 2 7 ’ 

(34) 

if I 7 I < y m (Eq. (i27l) i. as in the case of Fig. [4] (cos 2 9' is 
given by the inverse of (l34ll if 7 m < 7 < 7 r — 7 m ). Hence, 
in finite chains small yet finite side-limits together with 
a discontinuity will be exhibited by the concurrences Cy 
as h crosses h' s at a fixed orientation, as verified in the 
top right panel of Fig. [I] 

We can extend Eq. (1551) to general spin S > 1/2 by still 
considering the reduced states pfj arising from | 0 ^j), as 
those of two effective qubits, stemmin g fr om the single 
site states 1©^}, |0^), as done in ref, fill] for the TFF. 
The generalized expression is obtained replacing cos O' —> 
cos 2S 9' and sin 2 9' — > 1 — cos 4S 9' in (1551) . The negativity 
can be similarly evaluated D- 

The side-limits at h' s of the reduced state of m given 
spins, p ± (m), can be directly obtained from the exact 
side-limits (1551) of the full GS. They will be rank 2 mixed 
states (and not rank one states, i.e., pure states, as in 
h s ), of the form 




FIG. 4. (Color online) The concurrences Ci between spins i 
and i + 1 (top) in a FM (left) and AFM (right) finite spin 1/2 
XY chain with \ = Jy/Jx = 1/2, as a function of the mag¬ 
nitude h = |/i| of the non-transverse field at fixed orientation 
ra 7 = (siny,cosy) in the xz plane, with y = 0.0271 (FM) and 
0.3671 (AFM). For these values there is a single factorizing 
field |/i s | « 0.76 J x (Eq. (1191) ) in the FM case, determining a 
UGS, and two factorizing fields \h' 3 \ « 1.06| J x \ (Eq. (1261) 1 and 
|/i s | « 1.431 J x \ in the AFM case, corresponding to a NGS and 
UGS respectively. The insets depict the details in the vicinity 
of these fields, showing that all CVs vanish linearly at h a (Eq. 
m) and approach the finite /-independent side-limits (031) 
at h' s . All pairs are entangled in the vicinity of h a and h' s , 
remaining so between both fields in the AFM case considered. 
All labels dimensionless. 


P ± M = 

|9w)(ey|+|-e);)(-e%|±(|e%)(-e)G(-e"- m |e"- m )+fr.c.) 


2(1±(—@iv|©iv)) 


(35) 


where | ± 0 ) 0 ) denote the reduced states of the m spins 
in the Neel states | ± ©at) and (—©a/Qat) = cos"0', 
<-©5vl©Jv) = c°s m 0', with cos 2 9' given by (El for |y| < 
y m . The exact eigenvalues of p ± (m) are p ± (m) and 1 — 
p ± (m), with 


p ± {m) = 


(l + cos m 0 ')(l±cos n - m 0 ') 
2(1 ± cos" 9') 


(36) 


The spectrum of p{m) will then reduce to these two eigen¬ 
values for h —> h' s . These side-limits are independent 
of the choice of the m spins, i.e., the same for m con¬ 
tiguous or sep arated spins, as in the UGS of the trans¬ 
verse case [ll|. For general spin S , we should just replace 
cos0' —> cos 2S 9' in (1561) . 

The ensuing side-limits at h' s of the entanglement en¬ 
tropy S[p(m)] are then 

S[p ± ( m )] = ~V ± (m)\og 2 p ± (m)-[l-p ± (m)} log 2 [l-p ± (m)] 

(37) 

For sufficiently large m < n/ 2, the overlap (—0)010)0) 
vanishes and p ± (m) —> 1/2, S'[ / o ± (to)] —> 1. For m = 2 
we also obtain from (1551) the side-limits of the reduced 
state of a spin pair, which lead to the separation inde¬ 
pendent limits (1331) of the concurrence. 


C. Discussion 

In Figs. HED w e show illustrative exact results for a 
cyclic FM (left) and AFM (right) spin 1/2 chain of n = 12 


spins interacting through first neighbor XY couplings 
( J z = 0 ) with x = Jy/Jx = 1 / 2 , immersed in a non- 
transverse field, where all previous effects can be clearly 
appreciated and verified. The numerical results were 
obtained through diagonalization (notice that an exact 
analytic solution of the XY chain through the Jordan- 
Wigner fermionization [32) is feasible just for transverse 
fields 0 )- All quantities are depicted as a function of 
the scaled magnitude \h\/\ J x \ of the non-transverse field 
at fixed orientation in the xz plane (7 = 0.027T in the FM 
case, 7 = 0.367T in the AFM case). For these orientations 
there is a single NTFF h s in the FM case, determining 
a UGS, whereas in the AFM case there are two NTFF, 
the first one h' s corresponding to a NGS and the second 
one h s to a UGS. 

It is first verified in the top panels of Fig. [5] that while 
at weak fields just the first neighbor concurrence C\ is 
finite in the present FM and AFM cases, all concurrences 
Ci become non-zero in the proximity of the factorizing 
fields. As seen in the insets, in the vicinity of h s their 
behavior is correctly described by Eq. flUD, vanishing all 
linearly with \h — h s | for h — > h s ((3ij oc Sh±r] I 1 ^ (77 > 
1) in the case of Fig. [4] changing sign as h crosses h s ). 
On the other hand, in the AFM case they all approach 
the finite /-independent distinct side-limits (1551) at h' s 
(here cos0' ~ 0.92 and C~ ~ 0.11, C + ~ 0.049). Both 
factorizing fields appear successively as the field increases 
along orientations 9 < 7 < y m , leading to a rather broad 
interval of “long-range” pairwise entanglement located 
between h' s and h s , as appreciated in the right panel. 
It is also seen that all Ci exhibit jumps for \h\ < |/i' |, 
the last one at h' s , which reflect the n /2 translational 
parity transitions of the exact GS, as discussed below. 
We remark that while the side-limits (1551) at h' s diminish 
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FIG. 5. The block entropies S(m) = S[p{m)] of m contiguous 
spins (top) and the eigenvalues pi (entanglement spectrum) of 
the corresponding reduced states p(m) for m = 5, in the same 
FM (left) and AFM (right) systems of Fig. H At h s (UGS), 
S(m) and all but one (pi) of the eigenvalues of p(m) vanish, 
while for h —y h'f (NGS), S(m) approaches the finite side- 
limits (1371) (indicated for m = n/2), and two eigenvalues (pi 
and P 2 ) remain nonzero, with pi approaching the indicated 
side-limits (1361) . The insets depict again the details in the 
vicinity of the factorizing fields. 


as the number n of spins increases, the finite values of the 
Ci ’s in the vicinity of both h' s and h s persist for larger 
sizes. 

The block entanglement entropies of to contiguous 
spins are depicted in Fig. [5j Block entropies in XY or 
XYZ spin chains have been studied in detail just for 
zero or transverse fields JT§, [H, [33|, including also block 
Renyi entropies mm It is first verified that the von 
Neumann entropies S[p(m)] vanish essentially quadrati- 
cally in the vicinity of h s , whereas in the AFM case they 
approach the finite side-limits (EZD at the Neel factoriz¬ 
ing field h' s (here S'[p + (m)] ss 0.31 while S[p~(m )] = 1 
for m — n/2; notice from (13^1) that p~(m) = 1/2 for 
to = n/2 V n). In the FM case these entropies rapidly 
saturate as m increases for all non-zero fields, showing 
then a non-critical behavior, whereas in the AFM case, 
while above |/r s | they become small (< 0.01) and also 
rapidly saturate, below \h' s \ they are larger and show an 
appreciable dependence with block size. 

The behavior of these entropies can be better under¬ 
stood by means of the entanglement spectrum, shown in 
the bottom panels, where the eigenvalues pi of p(rn) for 
to = 5 are depicted. Results for other m > 1 are similar. 
In the FM case, there are three dominant eigenvalues, 
with pi close to 1, and the behavior of p -2 and p% re¬ 
sembles that of S(m): All eigenvalues except p\ vanish 
(quadratically) at h s . However, in the AFM case it is 
seen that for \h\ < |h'|, both pi and P 2 are significant 
and comparable, indicating roughly an approximate rank 
2 reduced state. When h — y h'f, p(m) becomes exactly 



h/J x h/\J x \ 

FIG. 6. (Color online) The scaled intensive magnetizations 
= JA 2 (S^)/n and m = |m| (top) and the scaled energy 
gap between the ground and the first two excited states (bot¬ 
tom) in the same FM (left) and AFM (right) systems of Fig. 
[4] Notice that m = 1 at the NTFF h s determining the UGS, 
as seen in the insets. The energy gap shows that the UGS is 
well separated from the first excited state, while the NGS is 
two-fold degenerate. All labels dimensionless. 


a rank 2 state and just pi and P 2 are nonzero, in agree¬ 
ment with Eqs. (l36l) (|35|) . The behavior is similar to that 
observed for a transverse field [l9| . As expected, the GS 
transitions taking place in this sector are clearly visible in 
both the block entropy and the entanglement spectrum. 
For \h\ > h' s there are just three dominant eigenvalues, 
with pi much larger than the rest, as in the FM case. All 
but pi vanish again at the second factorizing field h s . 

Let us remark that at zero field, the results for any 
entanglement measure in the FM and AFM XY chains 
with first neighbor couplings are strictly coincident, since 
the corresponding Hamiltonians can be transformed into 
each other by a local rotation of angle n around the z 
axis at all even sites, which does not affect entanglement 
measures. This fact explains the pronounced increase of 
the block entropies of the FM case as the field vanish, 
since they approach in this limit the higher AFM val¬ 
ues. This symmetry no longer holds for finite fields not 
pointing along the z axis. 

In spin 1/2 systems, the magnetization can be used as 
a separability witness: The quantity m = 2 JA KS^I/n 
satisfies to < 1 in any pure entangled state of such sys¬ 
tem, with to = 1 if and only if the pure state is com¬ 
pletely separable. For a state with TI, to = 2| M\/n, with 
M = (JA Si) the total magnetization. Hence to = 1 at 
the NTFF h g , as verified in the top panels of Fig. [5J 
entailing a non-monotonous behavior of to. for increas¬ 
ing fields, as seen in the insets. We have numerically 
checked that such non-monotonous behavior persists for 
larger sizes, indicating that it is not a finite size effect. 
Therefore, through a careful measurement of M or the 
associated susceptibility as a function of the applied field, 
one could be able to identify the NTFF h s . 

In the bottom panels of Fig. [5] it is seen that the UGS 
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FIG. 7. (Color online) The concurrences Ci between spins i 
and i + l in FM and AFM XYZ chains with J y /J x = 1/2, 
J z = 0.2\J X \ and J x > 0 (< 0) in the left (right) panel. The 
orientation of the applied magnetic field in each panel is the 
same as that of Figs. [4][6] The factorizing fields are now \h s | ~ 
0.52J3, in the FM case, where x = 0.375, and \h s \ ~ 1.39 J x , 
\h' s \ ~ 0.84TC in the AFM case, where x — 0.583. The side- 
limits at h' s and the linear vanishing of all concurrences 
at h s are again verified. 


at h s is non-degenerate and well separated from the first 
excited state, whereas the NGS at h' s is two-fold degen¬ 
erate. Actually, as seen from the energy gap and also 
from the magnetization and previous entanglement mea¬ 
sures, while no transitions are observed in the FM case, 
in the AFM case the exact GS exhibits n/2 transitions 
as \h\ increases at fixed 7 , the last one taking place at 
h' s . They correspond to “translational parity” transi¬ 
tions IGS^) —> IGS'^), with |G'5' =b ) the exact TI ground 
states, which satisfy TIGS' ± ) = ilGS^} (T is the one- 
site translation operator). These transitions are similar 
to those observed for transverse fields in both AFM or 
FM systems Games in, where they are related with 
spin parity transitions and also end at the corresponding 
transverse factorizing field m m, EH- Hence, h' s still 
represents, in the non-transverse case, a critical field for 
the finite system, indicating the passage to a different 
regime. 

Finally, we depict in Fig. [7] results for the pairwise 
concurrence in a chain with full XYZ couplings, a system 
which cannot be mapped to independent fermions even in 
the transverse case ( |32[ IS)- We have set J z = 0.21|, 
in both the FM ( J x > J y > 0) and AFM ( J x < J y < 0) 
cases, using the same previous field orientations. The 
behavior is quite similar to that of Fig. 01 with the GS 
translational parity transitions also present in the AFM 
case. One just notes the higher values of G; above the 
factorizing fields in both cases, and the closer side-limits 
at h' s in the AFM case (now C~ ~ 0.036, G + ~ 0.032), 
due to the different value of the anisotropy ratio y. The 
values at zero field are again still strictly coincident due 


to the same sign of J z . Results for the block entropy and 
entanglement spectrum for the finite case considered are 
also qualitatively similar to the previous results. 

IV. CONCLUSIONS 

We have first determined the general conditions for the 
existence of separable eigenstates with maximum spin at 
each site in general arrays with anisotropic XYZ cou¬ 
plings immersed in a non-transverse field. The set of 
factorizing fields can be characterized by the local fields 
orthogonal to the local alignment direction, plus arbi¬ 
trary fields parallel to the latter. We have next identified 
the possibility of a uniform non-degenerate separable GS 
in quite general systems of arbitrary spin, including FM 
and AFM-type chains and arrays, for fields parallel to a 
principal plane (Fig. [2]). The coupling range can be ar¬ 
bitrary, provided the anisotropy ratio \ is constant. In 
AFM XYZ chains with first neighbor couplings, this sep¬ 
arable solution coexists in field space with the Neel-type 
separable solution. 

We have also demonstrated, for arbitrary spin, that 
pairwise entanglement reaches full range in a finite array 
in the vicinity of the factorizing field determining the uni¬ 
form solution, with the concurrence vanishing linearly at 
this field. Full range is also reached at the Neel NTFF, al¬ 
though here it was shown that in finite cyclic even chains, 
the pairwise concurrence reaches finite side-limits in its 
vicinity, which were analytically evaluated. This NTFF 
was shown to correspond to the last parity transition of 
the GS in the finite cyclic chain. Block entropies were 
also analyzed and shown to vanish quadratically at the 
uniform NTFF, while reaching again finite (and analyti¬ 
cally determined) side-limits at the Neel NTFF in these 
finite chains. 

Present results and limits are also applicable to more 
complex systems, like dimerized chains and arrays EH EH, 
l36| . The recent possibility of performing quantum simu¬ 
lations of spin chains and lattices with tunable couplings 
through cold atoms in optical lattices f37l - f39l | or trapped 
ions |39HTll augments the potential of the present results. 
Such experiments could then provide valuable insights 
into the remarkable phenomenon of factorization and its 
relation with entanglement and criticality in finite many 
body systems. 
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